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ABSTRACT 

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid 
masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are 
relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management 
in a low-gravity environment. Specific issues include (i) the mechanics of liquid masses in large-amplitude 
motions, (ii) the influence of bounding surfaces on the motion and (iii) the ability of such surfaces to 
control liquid motion by wetting forces, especially when they are augmented by various surface treatments. 
Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, 
and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is 
divided into an ’outer 5 or inviscid potential-flow region and one or more inner, or viscous dominated, regions. 
Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow 
droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation 
uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior. 


INTRODUCTION 

From an engineering viewpoint, it is often necessary to maintain the location of a quantity of liquid 
in an enclosed vessel. This task is made simple in terrestrial gravity where the liquid in a partially-filled 
container can naturally be assumed to locate itself near the ’bottom’ of the vessel. In the absence of gravity, 
however, it is possible for a liquid mass to become disconnected and become intermixed with its own vapor 
in an a priori unpredictable manner. Similarly, when the stabilizing gravitational potential is not present, 
the center of mass of a liquid region may move a large distance when disturbed by relatively small applied 
forces, with influence on the overall dynamics of a space vehicle. 

Liquid motion in microgravity allows a more complete examination of surface tension and wetting forces. 
In normal gravity, capilllary forces are restricted to short length scales and their influence, while important, 
is often obscured by the concurrent high level of viscous damping. When terrestrial gravity is removed, the 
spatial extent of capillary-dominated regions is greatly magnified and more detailed measurement is possible. 

A strategy for control of liquid motion involves the location and surface treatment of container walls. 
The motion of a liquid front on a solid surface is governed in part by the energetics of the three interfaces 
that meet at the so-called contact line. A local force balance specifies the static contact angle between the 
phases there. Terrestrial experiments have indicated surprisingly high rates of damping of liquid motions 
which has been attributed to energy dissipation at moving contact lines [1]. For a surface whose energies 
vary with position, e.g. if the surface is ’contaminated, 5 it is possible for the contact line to be trapped in 
various local energy minima leading to a significant difference, or hysteresis, between advancing and receding 
contact angles [2]. Forced motion of the contact line on these surfaces will require a succession of small energy 
barriers to be overcome and may provide sufficient dissipation to prevent large motions. Quasi-static analyses 
by Schwartz & Garoff [3,4] considered various patterns of surface wettability. Some of their predictions were 
confirmed experimentally. Those results are only appropriate to very slow imposed motions. One goal of the 
present project is the extension to finite speed, where viscous forces must also be considered. Small scale 
roughness is believed to play a role similar to surface contamination [5]. 
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THE LARGE-SCALE INVISCID PROBLEM 


Large-scale oscillations of liquid masses involve, principally, interchange between the kinetic energy of 
motion and potential energy stored as elongation of the bounding surface. Solutions for wave motions of 
infinitesimal amplitude are well known. For example, the small-amplitude oscillation of a liquid droplet was 
predicted by Rayleigh in 1877 [6]. Large amplitude motions, on the other hand, are highly nonlinear and 
will, in general, require numerical solution. Sloshing motions of low viscosity liquids in tanks of moderate or 
large size may, away from the walls and contact lines, be considered to be irrotational potential flows. The 
potential <j> is found as the solution of Laplace’s equation 

V 2 <^ = 0, (2.1) 

subject to the kinematic boundary condition 

DF/Dt = 0 (2.2) 


on the free interface F(x, y, z, f) = 0. Here D/Dt is the substantial derivative. The solution must also satisfy 
the condition of no normal relative velocity at stationary or moving impermeable boundaries. Additionally, 
on the free surface, whose shape evolution needs to be found as part of the solution, we have the dynamic 
condition 

<j>t 4 (1/2 ){4> 2 x 4- <f>l 4 <t>l) 4 4 <r ^ 4 y = C(t) . (2.3) 

Here p is the constant liquid density, cr is surface tension, and C(t) is a function of time only. R\ and R 2 are 
principal radii of curvature. The applied surface pressure is ordinary taken equal to zero without loss 
of generally; here, however, it may be used to represent the far-field influence of the ’inner’ moving contact- 
line problem. Time-dependent solutions of (2.1) - (2.3) can be achieved using boundary-integral methods 
for two-dimensional or axisymmetric problems. These techniques find generalization in the so-called panel 
methods for three-dimensional problems. An example of a boundary-integral solution to a two-dimensional 
problem is shown in Fig. 1 where a limiting wave form traps a vapor bubble [7]. 

DROPLET MOTION ON SUBSTRATES OF MIXED WETTABILITY 


We treat the low-Reynolds-number motion of a thin Newtonian liquid layer onto a previously-dry surface. 
The mathematical model employs the long- wave or ’’lubrication” approximation. To leading order in the 
free-surface inclination, only a one-dimensional, unidirectional flow problem needs to be solved. Performing 
this integration across the thin dimension first results in a reduction of the dimensionality of the problem. 
Consistent with the lubrication hypothesis, many physical effects of interest can be modeled much more 
simply. Asymptotic derivations of the lubrication model, for two-dimensional problems with a free surface, 
are given by Benney [8] and Atherton & Homsy [9]. 


Within this approximation, integral mass conservation and the creeping-motion force balance are com- 
bined to yield 

(3.1) 


1 


h t = -V • Q = • (/» 3 Vp) 


p =-„v>*-n = -„v’ k - B (JL- _i_) (3.2) 

Here h is the coating thickness, p is pressure and p is viscosity. The V - operator is two-dimensional, and 
is written using the substrate coordinates, and t is time. The first term in (3.2) incorporates the thin-layer 
approximation to the surface, curvature while the second is a two-term model of disjoining pressure that 
allows contact-line motion, as described below. B and h * are positive constants. The form of the disjoining 
term depends on the material system. Teletzke et al give various choices including the one used here [10]. 
Several flow simulations, using a form of disjoining pressure, have appeared [10,11]. 

When a liquid advances onto a previously-dry substrate, the usual mathematical description of the 
motion fails in the immediate vicinity of the three-phase line, i.e. the line where the liquid, solid, and 
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the vapor meet. The standard boundary condition, when a viscous liquid is in contact with a solid, is no 
relative motion between the liquid and the solid, the so-called ’no-slip’ condition. It has been shown that 
solutions incorporating the no-slip condition lead to paradoxical results; the work required to move a liquid 
onto a dry surface is infinite [12,13]. Clearly the no-slip condition must be abandoned very near the moving 
contact point or line. Several slip models, involving a free parameter, have been proposed [14,15,16]. Other 
possibilities are a uniform precursor layer [17], a grid-slip model that allows motion onto dry regions, with 
the rate of advance dependent on the numerical mesh spacing [18], and shear-thinning rheology near the 
contact line [19]. Figures (2a) and (2b) show calculations for liquids advancing onto nominally-dry substrates. 
Figure (2a) shows the development of fingers or drip marks for gravity-driven flow on a vertical substrate 
[18]. Figure (2b), using a uniform precursor layer, simulates finger development on a spherical body. 

The inclusion of the the disjoining terms in (3.2) is another technique that will allow contact-line 
motion using a thin wetting layer. Consider the spreading of an initially axisymmetric droplet onto a 
nominally-dry substrate. The evolution equation is nondimensionalized by measuring h in units of ho, a 
typical film thickness. The substrate coordinates (x,y) are measured in units of L and the unit of time is 
T* = 3fiL 4 /(<rh 3 ). In dimensionless form, the evolution equation is 1 


ht = -v- (/» 3 vv 2 /») + cv • [(i - A) ™] 


(3.3) 


where 



(3.4) 


The now-dimensionless parameter h* « 1 and, for droplet spreading, when combined with the value of C\ 
prescribes the value of the equilibrium contact angle on a given substrate material. The energy per unit area 
of the substrate is composed of two parts. Within the lubrication approximation, the surface tension energy 
density is 

e (<7 > = io-V/i • Vh (3.5a) 

and the disjoining energy density is 


,(<*) 





(3.56) 


both in dimensional units, has a minimum at h = (3/4)h* and (3.5b) is normalized so that is zero 
there. It is possible to identify the value of B/h * 3 with the equilibirium contact angle. For a static situation, 
= e( d )(3h*/4) = 0 and = e^(h — ► oo) when h* is very small. A local force balance at the contact 
point gives 

<7COS 0 C = fsv - JSL = O - ( 36 ) 

which, consistent with the small-slope approximation, yields 




1/2 


(3.7) 


Here the subscripts S, L, and V have their usual meaning of solid, liquid, and vapor respectively. Since 
9 C is the material property to be reproduced, constancy of the quantity C/h * 3 is to be maintained as h* 
becomes small. Provided h* is taken sufficiently small, simulation results become independent of its actual 
value. Numerical tests reveal that a droplet will stabilize with a maximum inclination approximately given 
by (3.7), smaller values of h * producing closer agreement. In this sense the actual value of h * is not needed; 
the only material property required is the apparent static contact angle. 

Few surfaces are ever ideally clean. Moreover, one may wish to purposely place patterns of contamination 
on a substrate in order to arrest liquid motion. Equation (3.3) is capable of simulating such a situation by 
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prescribing different values of C in particular substrate locations. Figure 3 shows the result of a time- 
dependent simulation where a drop of non-equilibrium shape is placed approximately at the center of a 
’’cross” of poor-wetting material. The contact angle, using (3.7) is twice as large on the cross, compared to the 
surrounding field. Reflection symmetry conditions are used at the boundaries of the computational domain. 
Figure 4 contains contour plots of the same simulation with the wettability boundaries superimposed. Rather 
realistic features are exhibited by the simulation; apparent advancing and receding contact angles are speed- 
dependent and the motion proceeds by alternating slow and rapid events. This is readily seen in Figure 5 
where the variation of the integrated energy components, using (3.5), is shown. This is similar to the result 
of the quasi-static analyses of Schwartz k GarofF [3,4] who termed the process ” sticking, stretching, and 
jumping” in their model of contact-angle hysteresis. The simulation ends with the droplet finding a stable 
energy minimum by merging with its neighbors in one corner. The rate of decrease of the total energy in 
Figure 5 can be shown to be equal to the rate of viscous dissipation. Notice than when the motion has 
ceased, the energy is approximately equally-partitioned between its two components. 

These computations use an alternating-direction-implicit (ADI) algorithm based on the scheme of Peace- 
man k Rachford [20], generalized to the present higher-order diffusive system. Nonlinearity is treated simply 
by evaluating prefactors at the old time level. The algorithm is quite stable and time stepping is adaptive, 
with larger steps being used when the surface shape is changing slowly. Steps can be as large as a factor 
of 10 7 bigger than the limiting step size for an explicit method. Basic features of the simulation have been 
verified to be invariant under temporal and spatial refinement and the choice of values for h * . Additional 
tests of convergence are currently in progress. The run shown here, using a 100 x 100 finite-difference grid, 
took about 30 minutes on a Pentium 100 PC configured as a UNIX workstation. 

We intend to consider other disjoining pressure models to assess the sensitivity of flow behavior to 
this choice. Physical experiments on surfaces of controlled wettability can be compared with numerical 
predictions. For relatively large drops, in the laboratory, gravity effects will need to be considered. These 
can be inserted in the model without difficulty. 
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FIGURES 



Fig. 1: A capillary-dominant progressive Fig 2; Flow simulations using two different wetting models, (a) A liquid 

free-surface wave calculated by a boundary- mound on a vertical wall shows characteristic fingering patterns [18]. 

integral method. A bubble is trapped [7]. (b) Drainage on a sphere, using a uniform precursor layer. 




Fig 3: Simulation of droplet motion on a mixed-wettable surface. Parameter values are Cl = 0.0002, C 2 = (1/4) Ci, 
h* = 0.02 . The grid size is 0.01 L where L is the width of the window. Times are t/T* = .001, .5, 4., 5. 6., 7. The 
initial profile was h = (3/4) [ h* + exp (-64 r^} ]. 
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Fig 4: Contour plots for the profiles in Fig. 3 . The wettability boundaries are also shown. 



Fig. 5: Variation in integrated energies with t/T*. The motion proceeds by alternating fast and slow events. Energies are 
measured in units of a h 0 2 . 
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